function [] = number_system(n,r,L,M,x)
    disp('----------------------------------------------------------------');
    fprintf('** System no. %d\n',n);
    disp(' ');
    lambda = 1; % positive radix system

    disp('* Larges representable number');
    largest_representable_number(r,L,M,x,lambda);

    disp('* Number of representable numbers');
    x_len = length(x);
    fprintf('Number of digits: %d => %d^(L+M) = %d^(%d+%d) = %d\n', x_len,x_len,x_len,L,M,x_len^(L+M));
    disp(' ');

    disp('* Unit of last position, ulp');
    fprintf('r^-L=%d^-%d=%.32f\n',r,L,r^-L);

    disp(' ');
end

